Iterated dynamical maps in an ion trap. 
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Iterated dynamical maps offer an ideal setting to investigate quantum dynamical bifurcations and 
are well adapted to few-qubit quantum computer realisations. We show that a single trapped ion, 
subject to periodic impulsive forces, exhibits a rich structure of dynamical bifurcations derived from 
the Jahn- Teller Hamiltonian flow model. We show that the entanglement between the oscillator and 
electronic degrees of freedom reflects the underlying dynamical bifurcation in a Floquet eigenstate. 



Iterated area-preserving maps are better adapted to 
physical implementations of quantum information pro- 
cessing than Hamiltonian flows. In fact one of the 
first quantum algorithms, the Grover search algorithm [H, 
may be regarded as a quantum description of an area- 
preserving map. Typically, a quantum algorithm consists 
of a sequence of elementary unitary operations on one or 
more two-level systems. In the Grover search algorithm a 
simple product of unitary operators is iterated. While it 
is possible to simulate an arbitrary Hamiltonian flow as 
a quantum circuit, area-preserving maps are simulated 
more directly as an iterated sequence of unitary gates. 
Originally introduced by Poincare, iterated maps have 
formed the core of studies in quantum chaos for many 
decades [3]. The model reported here is a simple example 
of a non trivial area-preserving map, based on the Jahn- 
Teller model d and is naturally adapted to an ion trap 
implementation of quantum information processing. 

The Jahn- Teller model describes a class of systems in 
which one or more particle coordinates are coupled to 
a two level system^. The E ® e Jahn- Teller model is 
a minimal description in which two harmonic oscillator 
coordinates are coupled to a two level system. We may 
write the Jahn -Teller E ® t Hamiltonian as 

Hn = AS:, + — {pI+pI)-\ ^ill+ll) + ^^1xSx+XyqySy 

. (1) 

where {q^, <ly,PxiPy} a-re the conjugate position and mo- 
mentum operators in the harmonic oscillator space, and 
{sx, Sy, Sz} are su{2) operators acting on the states of the 
two-level internal degree of freedom. 

The classical description of this model is a Hamilto- 
nian flow in the phase-space of the Cartesian product of 
the phase plane of the oscillators and the spherical 
phase space of the two level system 5^, known as the 
Bloch sphere. The quantum model exhibits the conical 
intersection that has assumed importance in various bio- 
physical models of visionfs^. It is well known that the 
E ^ e Jahn- Teller model displays a classical bifurcation 
of fixed points as the coupling strength is varied and that 
the quantum analog displays a maximum in the entangle- 
ment at the same value of the coupling strength 6, 7] as 
the classical fixed point bifurcation. At this point there 
is a morphological change in the nature of the ground 



state. This is the few body analogue of a quantum phase 
transition in a true many body system. 

We do not consider the Hamiltonian flow model, but 
rather a related model, which we will call a kicked E ® e 
model that is better adapted to an ion trap realisation. 
Instead of a time-independent Hamiltonian, we will con- 
sider a strongly time dependent version of the Hamilto- 
nian. We assume that 
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A hamiltonian which is periodically varying in time is 
most naturally described in terms of the Floquet opera- 
tor, a unitary operator that maps the dynamics over one 
period. 

Two kicking interactions are periodically applied in- 
stantaneously and consecutively, after which there is a 
period of evolution r under the Hamiltonian flow associ- 
ated with the confining potential. 

The unitary operator which discretely maps states over 
each period is the Floquet operator (setting h= 1) 

U = eiip{—iHoT) exp{—iXHx) exp{—iXHy), 
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Ho = As, + ^(pI +pI) + ^{ql + ql), 



(2) 



The unitary map ^ is the kicked equivalent of the E ® 
e Jahn-Teller model. We will fix the frequency of the 
harmonic potential, lu, the internal splitting energy. A, 
the periodicity of the kick r and adopt units such that 
TO — uj^^. This leaves the kick coupling A as the tunable 
parameter. 

For each system observable A evolving under a Hamil- 
tonian H the time evolution is d{A)/dt = i{[H,A]). 
From this expression, we can derive operator differen- 
tial equations for the evolution of all seven operators un- 
der each of Hq, H^, Hy. According to Ehrenfest's theo- 
rem [8] , in the classical limit expectation values of oper- 
ators approach classical variables. Hence the differential 
operator equations become a set of classical equations 
in this limit. Integrating these equations over one pe- 
riod T gives three maps which we then compose into a 
single area-preserving map on the classical phase space 
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M(M'* X 5^) R4 X 52 where 5^ denotes the Bloch 
sphere. 

The operator time evolution equations for 
<lx, Px, Py, Sx, Sy, Sz can be determined under 
each of exp(— jiJo), exp(— iiJ^;), e:xjp{~iHy). Taking 



these to be classical differential equations, in each case 
we can integrate these to obtain three maps for the 
corresponding classical variables. Composing these three 
maps yields 



Px 
Sx 



XSy 



{px — A{si: cos(A(jy) + sin(Agj;)}) sintj + q^ cosu, qy 
(px — )^{sx cos{Xqy) + Sz sm{\qy)}) cosLo ~ qx sinuj, Py {py 
(cos A cos(qyA) — sin A sin(qa; A) sin(gyA))sa: — sin A cos(qa;A)si, 

+ (cos A sin(gyA) + sin A sin(ga:A) cos((j'yA))sz 
(cos A sm{qxX) sin(qyA) + sin A cos(gyA))sa: 

+ cos A cos{qxX)sy + (— cos A sin(ga:Aa) cos(gyA) + sin A sin(gyA))s 
— sin((jyA) cos(ga:A)sa; + sm{qxX)sy + cos(gyA) cos{qxX)sz 
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Py — XSy) sintj + qy ( 
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(3) 
(4) 

(5) 

(6) 
(7) 



Our first step is to solve for fixed points of the map by 
calculating solutions to M{x*) = x* . The fixed points 
correspond to periodic orbits of the system in phase 
space. There are trivial fixed points at the equilibrium 
of the harmonic oscillator where the pseudo-spin is ei- 
ther aligned or anti-aligned with the z-axis. These fixed 
points are the only ones present when A = 0, one of which 
is stable, 

Sz = -1/2, Sx = Sy = Qx = Qy = Px = Py = 0. (8) 

while the other (s^ — 1/2, all other co-ordinates equal to 
zero) is unstable. As the kicking coupling A is increased 
a bifurcation of the stable fixed point occurs when 
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8tan(tj/2) 
cot(A/2) ± 1 



(9) 



where uj = wr, A = Ar, are dimensionless variables. De- 
pending on the values of uj, A, the above equation can 
admit either 0,1 or 2 solutions. The smallest value of 
A for which a solution exists corresponds to a pitchfork 
bifurcation where the fixed point ([5]) becomes a saddle 
point and two new stable fixed points emerge, as illus- 
trated in Fig. [H When there is a second solution to ([9|) 
a second pichfork bifurcation occurs at the origin, with 
the saddle point becoming a local maximum and two new 
saddle points emerging. 

For A = we can associate the stable fixed point in 
phase space with the ground state of the quantum sys- 
tem. As A is varied we can (numerically) track this eigen- 
state of the Floquet operator ([2]), which we will refer to 
as the Pseudo Ground State (PGS), denoted \ipg),- Anal- 
ogously we define the first Pseudo Excited State (PES), 
denoted IV'e)- Below we demonstrate how the bifurca- 
tions in the phase space signal distinctive qualitative 
changes in the PGS, providing insights into a one-body 
analogue of quantum phase transitions which occur in 
many-body systems. 



For the oscillator space we choose the basis states 
which are simultaneous eigenstates of the number op- 
erator N — Px+pl+ll + Qy and the angular mo- 
mentum operator Lz — qxPy — QyPx, which we label as 
\N,l), N >l> -N. The basis |-|-), |-) for the internal 
two-level system are eigenstates of Sz with eigenvalues 
±1/2. The unitary map ([2]) is invariant under the parity 
transformation 11 — —iexp{iTTJz) where Jz = Lz + Sz is 
the total angular momentum operator. Since 11^ = / the 
Hilbert space splits into two parity classes given by 



O-{|2fc,0|- 
i? = {|2fc,0h 



\2k- 
\2k- 



1,01- 



k e N}, 
k e N}. 



As the state space of the system is infinite-dimensional, 
numerical diagonalisation requires a choice of truncation 
Nt- We take Nt = 18 giving the dimension of the Hilbert 
space to be 342, which is expected to give reliable results 
in the weak coupling limit since for A = the ground 
state is |0, 0) |— ) e O. As A is incrementally increased, we 
numerically solve for the PGS (in the subspace O) with 
an adaptive procedure which requires that the increment 
i5A is chosen such that 

\{^g{x)\ijg{x + sx))\<om. 

Husimi functions _9,] allow for the representation of 
quantum states in Hilbert space as a density in a clas- 
sical phase space, thus providing means for comparison 
between classical and quantum systems. They are de- 
fined in terms of coherent states 
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(10) 

where = (q^ + ip^)/^/2 and \n^) are number eigen- 
states of — p^ + q'^. The explicit form of the Husimi 
function is 

Il{tlj,ax,ay) = Tr{TTs{\tl)) {ip\)\ax, ay) {ux, ay\) 
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FIG. 1: Fixed point bifurcation of the classical system. 

Phase space cross sections in the — q.y plane are shown for 
the parameter values uj = 7r/60, A = 2 acot(2), giving the first 
bifurcation coupling as A^i = 0.26 and the second bifurcation 
at the origin when \b2 = 0.46. (a) A = 0.15 < Ati: There 
is a single stable fixed point located at the origin. Here the 
fixed point structure is invariant under any refiection axis 
which passes through the origin, (b) Af,i < A = 0.32 < Ai,2: 
Three fixed points occur with a saddle point at the origin 
and two other fixed points which are stable. There are only 
two refiection axes passing through the origin which leave the 
fixed points invariant. At A = A62 = 0.46 the fixed point at 
the origin changes from a saddle point to a local maximum, 
and two new saddle points emerge. 




(a) (b) 

FIG. 2: Husimi function cross sections, u — tv/GO and 
A — 2acot (2): (a) The one-dimensional cross section of the 
Husimi function is shown along the section = qy, Px = 
— tan(aj/2)g^ as the coupling A is increased. For small values 
of A the Husimi function is localised near q^ = 0. For A > 
Aft — 0.26 the Husimi function becomes first delocalise and 
again localised in two regions of the g^^-axis, indicating the 
wavefunction is a linear combination of localised states which 
we associate with the fixed points in phase space, (b) Cross 
section of the Husimi function in the qx — qy plane is shown for 
the state \tp+) = (IV^s) + |^e))/V2 when A = 0.32. The Husimi 
function is localised at co-ordinates close to the value of the 
fixed point in phase space, supporting the approximations 
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which exists in the four-dimensional oscillator phase 
space. Above, Tr^ is the trace over the two-level sub- 
system. 

In Fig. [2] the position space cross section of the Husimi 
function for the PGS is shown as a function of the cou- 
pling A. A signature of the pitchfork bifurcation of the 
trivial fixed point is apparent. A state that is highly 
localised at the origin of the oscillator space splits into 
branches associated with each of the classical stable fixed 
points after the bifurcation. For A < At, we expect the 
wavefunction of the PGS to be localised at the stable 
fixed point in the classical picture; i.e. the PGS is ap- 
proximately 

IV's)«|o,o)|-) 

where |0,0) corresponds to the harmonic oscillator 
ground state and |— ) is the eigenstate of Sz with eigen- 
value — 1/2. After the first bifurcation we expect the 
PGS to be in a linear combination of two states localised 
at each of the two new stable fixed points in the classical 
picture. A parity invariance of the Floquet operator, de- 
termines the phase relationship between these two states 
with the result that the PGS is approximately 

IV-s) ~ {\a^,ay)\n) - \ - a^,, -ay)\h')) (11) 

and the first PES is approximately 

IV'e) « (|a^,ay)|n) + \ - a.^,-ay)\n')) . (12) 

The state |n) = cos(e'/2)|+) + e'"^ sin(6'/2)|-) is deter- 
mined by the pseudo-spin component n — sin(0) cos(0)i-|- 



sin(6') sin(0)j 4- cos((?)k of the classical fixed point, and 
n' is the parity transformed pseudo-spin component ob- 
tained from (j) (j) + TT. 

By taking even and odd combinations of the states 
(|llll2p . we can test this approximation for A > A;,. It is 
expected the two combinations give states localised near 
each of the classical fixed points respectively, (although 
with some discrepancy due to mixing of the oscillator 
reduced density matrix which is absent in the approxi- 
mation Ea. (|12p . The Husimi function of the even com- 
bination (IV'g) + |V'e))/'\/2 of the numerically determined 
PGS and first PES is shown in Fig. [U^b), confirming that 
localisation occurs. 

As the PGS is a pure state the entanglement between 
two subsystems A and B can be quantified by the von 
Neumann entropy of the reduced density matrix 

S{pa) = -TripAlogipA)) 

where pA = TrA{p)- In this manner the entanglement 
between the two-level system and the oscillator degrees 
of freedom, or between an oscillator degree of freedom 
and the rest of the system, can be calculated from the 
numerically determined PGS. To calculate the entangle- 
ment between the two oscillator degrees of freedom we 
use the reduced density matrix po = Trs(\^){'ip\). which 
is generally a mixed state. In this instance the entangle- 
ment may be quantified through the logarithmic negativ- 
ity fl^ 

En{po) =logN II P^-" 111 

where the norm is the trace norm II A ||i=TrVAFAand 
p'^^ denotes the partial transpose over subsystem A. 
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FIG. 3: Entanglement measures, uj — n/60, and A = 
2acot(2): The first derivatives of three entanglement mea- 
sures for the PGS are shown as a function of A. From top to 
bottom the measures are the von Neumann entropy between 
the two-level and oscillator subsystems, the von Neumann 
entropy between the subsystem of a single oscillator degree 
of freedom and a pseudo-spin/oscillator pair subsystem, and 
the logarithmic negativity of the mixed-state reduced density 
matrix of the oscillator degrees of freedom. In all cases the 
derivative of the entanglement measure has dominant support 
over the quantum crossover regime Ai,i < A < A(,2. 



The model of this paper might be realised in a cylin- 
drical Paul trap. In order to ensure that the motion 
in the x — y plane is well confined we need to ensure 
that the secular frequencies in this plane smaller 
than the frequency in the z— direction. Using the stabil- 
ity diagram for this kind of trap (see figure 1 in [ll'l) we 
can operate at trap parameters such that ojx^y 
without significant micromotion in the x — y motion. Al- 
though this requires operation close to the stability edge 
it should be achievable. In order to couple different com- 
ponents of the Pauli matrix to the vibrational degree 
of freedom we can make use of the Raman scheme for 



state dependent displacements introduced by Monroe et 
al.[l3|. Consider first that the Raman beams are directed 
along the a;— axis of the trap but the same arguments hold 
for the y— axis Raman pulses. The effective interaction 
Hamiltonian for the Raman pulses is — x<lx<^z where 
CTz = |e)(e| — \g){g\ with |e), \g) are the ground and ex- 
cited states of the relevant two-level electronic transition. 
In other words we have a pseudo-spin realisation with 
1+) ~ |e), |— ) — \g). If we 'sandwich' this Raman pulse 
by two laser pulses tuned to the carrier transition with an 
appropriate phase choice, the displacement can be made 
to depend on any component of the Pauli matrices. For 
example the pulse sequence. 



(13) 



will achieve the desired result for the x— kick. 

We are now in a position to consider the experimental 
signatures of our results. As a specific example, we will 
focus on detecting the entanglement implicit in the state 
just beyond the first bifurcation, approximately given by 
in Eq. (fT2|) . This can be reached adiabatically starting 
in the state in |0, 0,)|— )and using a sequence of pulses 
with gradually increasing coupling strength. At the end 
of a sequence the probability to detect the atom in the 
excited state is simply 



P(+) = cos2(6l/2) 1-e 



(14) 



where 9,ax,0iy are determined by the classical fixed 
points. Sampling this distribution thus indicates the sup- 
port on the classical fixed points. Readout of the excited 
state is easily done in ion trap realisations using a cy- 
cling transition ll| . The experimental verification of the 
entanglement that results after the bifurcation is there- 
fore quite achievable with current technology. 
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